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Abstract 



This work investigates whether the observed 1 h exotic 
mesons are molecular states. We first use a potential model to 
calculate the spectra and lifetimes of the / (980) and a (980), 
taken to be loosely bound molecular states of KK, then 
apply the same scenario to the 1 h exotic states 7Ti(1400) 
and 7Ti(1600), assuming them to be 7r?y(1295) and 7r?7(1440) 
molecules respectively. We derive the effective potential in 
the framework of field theory at the hadronic level. Our re- 
sults indicate that the present data on 7Ti(1400) and 7Ti(1600) 
rule out the specific molecular ansatz. We show that the life- 
time of a loosely bound heavy-light molecule with enough 
angular momentum is fully determined by the lifetimes of its 
constituent mesons. 
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1 Introduction 



New gluonic matter, such as glueballs and hybrids, are predicted by the general theory 
of QCD. But after more than three decades of effort, gluonic excitations have not been 
discovered unambiguously. The discovery of the J PC = 1 h isovector resonances 7Ti(1400) 
and 7Ti(1600) is great progress along these lines, because it may imply that new gluonic 
matter has eventually been observed. The E852 Collaboration at BNL observed 7Ti(1400) 
in r/Ti in the reaction ir~p — > r/n^p, with a mass and width of 1370 ± 16^|q MeV and 
385 ±40i^ 5 MeV respectively 0. Independently, 7^(1400) has been observed in r/ir with 
M = 1400 ± 20 ± 20 MeV and T = 310 ± 50±|g MeV [|], or M = 1360 ± 25 ± 20 MeV and 
r = 220 ±90 MeV |3j by the Crystal Barrel Collaboration. The other resonance 71^ (1600) 
has been observed by E852 in pir in 7t~p — > 7t + 7t~tt~p with mass 1593 ± 8±|f MeV and 
a width of 168 ± 20i^° MeV [|; in rj'ir in n'p rf^-p with M = 1597 ± lOi^o M eV 
and T = 340 ± 40 ± 50 MeV ||; and most recently in a preliminary analysis in bin in 
Ti- p un-7i p with M = 1582 ± 10 ± 20 MeV and T = 289 ± 16 ± 27 @. The 7^(1600) 
is observed decaying to the same final states by VES 0. 

A 1~ + meson cannot be a meson (qq), a gluinoball (gg) or a squarkball (qq), because it 
is J PC exotic H. The observed 1 h mesons therefore should be considered as the hybrid 
or glueball candidates which are predicted by QCD. Since glueballs are isoscalar this 
possibility is excluded. The 7r!(1600) is consistent with expectations for a hybrid meson |], 
|TU|j , while the 71"! (1400) is not || [Dj]. Discovery of a hybrid is exciting, but identifying 
the observed 1~ + as the long-expected gluon excited states needs careful further study. 
Possible alternative hypotheses about their structure have been discussed [Kj. There was 
an attempt to explain the 71"! (1400) as a coupled channel effect dynamically generated 
in 7ii] and ttt]' scattering |y|. Donnachie and one of us (Page) suggested that the peak 
at 1400 MeV is a consequence of the interference of a 1600 MeV resonance with a non- 



resonant Deck-type background with appropriate phase |L3|. So far, this scenario has 
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only been argued to account for the E852 data, and has not been extended to the Crystal 
Barrel data. Also, aspects of the proposed mechanism have been challenged on theoretical 
grounds JTT], IT3]. However, it has been argued that there is no need for more than one 



resonance in the 1.3-1.6 GeV mass region when rjir, 77 tt and p-n data are considered JT5 



Before one can conclude that the observed 1 h resonances are gluon excited states, 
even more possibilities should be explored. Among the possible alternatives, we explore 
the hypothesis of molecular structure. For 1 h quantum numbers there may be many plau- 
sible meson-meson bound states and their mixtures. Since the binding energy for molecu- 
lar states is not very large, the sum of the masses of the constituent mesons should be close 
to the mass of the exotic state. According to the Particle Data Group jTjJ , one of the favor- 
able combinations are that 7Ti(1400) is composed of tt (0 h ) and 7/(1295) (0 h ) with total 
mass 771^ + m ?? ( 12 95) « 1430 ~ 1438 MeV, and 7i"i(1600) is composed of tt and 7/(1440) with 
total mass m n + m r) ( 1440 ) « 1535 ~ 1620 MeV. The constituent mesons in these molecules 
have one unit of relative angular momentum ("P-wave"). Since IV1295) = 53 ± 6 MeV 
and T ?? (i44o) = 50 ~ 80 MeV are not too wide, the 77(1295) and 7/(1440) build observable 
bound states with the pion, i.e. they can be constituents of molecular states with total 



widths of a few hundred MeV. The possibility has also been discussed elsewhere |L7j . In 
this work, we investigate whether the present data on 7Ti(1400) and 7i"i(1600) support this 
specific molecular hypothesis. Other possibilities for 7i"i(1400), not discussed in this work, 
are a b\ and tt with no angular momentum between them ("S-wave"), denoted (61^)5. 
Other S- and P- wave molecular possibilities are (fiir)s and (K*K)p. For 7i"i(1600) the 
possibilities are (/ 1 (1420)7r) s , (/i(1510)7r) 5 , (puj) P and (p(1450)7r) P . 

We use a potential model to calculate the spectra and lifetimes of ?Ti(1400) and 
7i"i(1600), and then compare them with data to determine whether the molecular hypoth- 
esis is reasonable. We derive the potential between the constituent mesons by writing 
down the elastic scattering amplitude of the mesons in momentum space and then car- 
rying out a Fourier transformation with respect to the 3-momentum of the exchanged 
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mesons to obtain the potential in configuration space []I8]. Substituting the potential 
into the Schrodinger equation, we can obtain the binding energy and wavefunction at the 
origin (or its slope) of the bound molecular state. The effective couplings are obtained 
from data by assuming SU(3) flavour symmetry. To account for off-shell effects in the 
calculation of the amplitude, certain form factors are introduced. We write down the 
inelastic scattering amplitude of the mesons in momentum space, from which the lifetime 
is calculated. 

To test validity of the mechanism employed, we first apply it to study /o(980) and 
do (980), which, in some interpretations, are interpreted as molecular states The 
molecular possibility has some support from lattice QCD [p0[| . As discussed above, we 
derive the potential between K and K, calculate the spectra of / (980) and a (980), 
obtain all necessary parameters and then calculate their lifetimes. The results are in good 
agreement with data. Thus we are convinced that the approach we adopt is applicable to 

1 h if it is a molecular state. 

The paper is organized as follows. After this introduction, we present most of the 
formulation and some technical details in Section 2. In Section 3, we report our numerical 
results and the last section is devoted to our discussion and conclusion. 

2 Formulation 

2.1 Exchange mesons 

Consideration of composites of two quarks and two antiquarks is complicated by the 
fact that different colour neutral combinations are possible. The first possibility is that 
no subcomponents of the system is colour neutral, commonly referred to as a four-quark 
state. The second possibility is that one quark-antiquark subcomponent is colour neutral, 
forming a meson, and that the same happens to the remaining subcomponent. This 
composite of two mesons is referred to as a molecule. The latter possibility will be the 
subject of this work. In the limit of a large number of colours N c in QCD there are no 



four-quark states, but only meson-meson molecules pT| . The molecular picture becomes 
more accurate if the mesons contained in the molecule are separated far from each other 
relative to the size of the constituent mesons. This would happen when the binding energy 
of the two mesons is small relative to the quark-antiquark energy within each individual 
meson (which approximately equals the mass of the meson when the quarks are light). The 
theoretical calculation of such a small binding energy of mesons is complicated because 
the binding energy is the result of delicate cancellations between attractive and repulsive 
interactions, each of which individually corresponds to much larger energies than the total 
binding energy. This follows because individual interactions are governed by the QCD 
scale Aqcd, which is much larger than the binding energy. The calculation of the binding 
energy has been done in lattice QCD |20| , |22f , QCD sum rules ||23|| , quark-antiquark pair 
exchange models |L4|, |4j], four-quark potential [19|, |25|, flux-tube |26| and other p7fl 



models, and molecular models with meson exchange [2S|. We adopt the latter. Meson 
exchange models should be most accurate when the binding energy of the constituents 
is low. This is because in this case the constituents are typically separated far apart, so 
that the cost that confinement imposes on exchanging coloured objects at large distances 
means that the exchange is more likely to be a colour singlet, i.e. a meson. Consideration 
of a system with low binding energy, i.e. of two mesons separated far apart, yields the 
result that exchange mesons with a long range would be most important for the calculation 
of the binding energy. Long range exchange mesons correspond to low mass m q exchange 
mesons, taken to be either the ir, o (also called the /o(400 — 1200)), i], p, lu, rj', / (980), 
ao(980) or in this work. An exchange meson can be neglected when its range ~ — is 



small with respect to the r.m.s. radius r rms of the molecular state, i.e. when — - — is 

771 q Trms 

small. 

The general formalism for meson exchange is reviewed in Appendix [A[ For the re- 
mainder of this section we describe the theoretical formulation used in our calculations. 
Because the momentum of the exchanged meson in t-channel is space-like while it is 
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time-like in s-channel, the propagator is expressed as 

1 



, t — channel (1) 

q 2 — m 2 



s — channel, (2) 



q 2 - m 2 q + iT m m q 

where the Breit-Wigner form is taken and T m is the total width of the exchanged meson. 
At vertices, we use the standard strong coupling forms JHJ) and the proper Clebsch- 
Gordon coefficients to manifest the isospin structure of the vertices. The couplings take 
the values extracted from corresponding decay widths. There are three types of decay 
modes relevant to this work: 

0+ ->0 - (T (S- wave), 1" -> (T(T (P - wave), 1~ -> 0"1" (P - wave). 

The decay widths are 

; ' W\\ (3) 



(2J+ 1) 8vrm 2 

where I is the magnitude of the 3-momentum of the decay product in the CM frame, m 
and J are the mass and internal angular momentum of the decaying meson, and M is the 
amplitude, with 

(i) \M\ 2 = gl for 0+-^ _ _ , 

{ii) \M\ 2 = gl{ Pl - p 2 r( Pl - p 2 y £ e x ;el for 1" - 0~0-, 

A 

(iii) \M\ 2 = gie a ^e al ^^P a P a '( Pl - p 2 y( Pl - p 2 y' 



xE^i<*E4<* fori" -0-1- (4) 

Ai A2 

where is the corresponding strong coupling, P is the 4-momentum of the decaying 
meson, pi,p 2 are the momenta of the outgoing mesons, and eA is the polarization vector 
of the vector meson. We note that all the mesons which are connected to the vertex 
are on their mass shells. Using Eqs. |], the coupling constant gi is calculated from the 
experimental width. If the experimental width is unavailable, <?j is obtained by assuming 
SU(3) flavour symmetry for the coupling. 
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Since the propagators are off-shell, we introduce form factors for the effective couplings 
g,i at the vertices. The usual form factors are expressed as |28|, |30] 



A 2 - m 2 



(5) 



A 2 - q 2 ' 

at each vertex with A is an adjustable constant which models the off-shell effects at the 
vertices due to the internal structure of the mesons. 



2.2 Effective potential between two constituent mesons 



Following the standard procedure given in ref. [0 and using the effective vertices, we can 
write down the scattering amplitude between the two constituent mesons in momentum 
space and convert it into configuration space by a Fourier transformation. This procedure 
was used to derive the effective potential between a quark and an antiquark [jy|]. We 
consider all pseudoscalar, scalar and vector meson exchanges in t— , s- and u-channel, and 
indicate explicit expressions for the processes that need to be included in our calculations. 
In the center of mass (CM) frame of the bound state, with 
(i) + the exchanged meson in t-channel, 



91 



'A 2 — m 



2\ 2 



1 



q 2 + m 2 



A 2 + q 2 I 4mim2 



1 - 



2mj 



2ml 



(2p-q-q 2 )( T ^+ ' 



Ami 



Am? 



(ii) 1 the exchanged meson in t-channel, 



9 2 2 



'A 2 



m 



2\ 2 



q 2 + m 2 y A 2 + q 2 



p ■ q 



q z 



m\m2 m\m2 Am\m<i' 



(6)' 



(7) 



(Hi) + the exchanged meson in s-channel, 



<7?(A 2 



m 



2\2 



Am±m2 



1 



/«2 



hp 



,12 



m?)(a? - A 2 ) 2 



m 2 ) 2 (a 2 



A 2 ) 2 



26p 2 



(6 2 - 26)p 2 



(a' 2 - m 2 )(a 2 - A 2 ) 3 2a 2 (a /2 - m 2 )(a 2 - A 2 ) 2 
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where we introduce the variables 



a 2 = 


(mi + m 2 ) 2 


b = 


(mi + m 2 ) 2 


m 1 m 2 


a' 2 = 


T 2 m 2 
a 2 + \* m * 
cr — m 


b' = 







(9) 
(10) 

(11) 
(12) 



where q is the 3-momentum that is exchanged. The momenta of constituent mesons with 
mass mi or m 2 in the initial state is respectively denoted by p and —p. After the Fourier 
transformation, we obtain the potential forms in configuration space: 



[i) + in t-channel, 
V(r) = 



9l 



167rmi7«2 



+ 



2\ a -Ar 



mj + m 2 . e mqT - e Ar 



2A 

,2 , _2 1^.2 



2m 2 m 2 



~2 + m l+ m 2 { m q ~ 



2\„-Ar 



ig 2 {A 2 -m 2 Q ) 2 



1m\m\ 2A 

— U(m ? r + l)e" m « r + 5(Ar + l) e - Ar 

2mi?7i2 

+(C-5A 2 ) e - A Vl fJL + 7 l ^ p r 

J \Am\ Am 

3 2 (A 2 -m 2 ) 5 



4m 2 , / 47rr 3 



4mim 2 



-Ar 



— — + Am 2 + 5m 2 - 

A 9 r 9 

-Ar 



+(C-BA 2 )m 2 



q A 



1 1 

+ 



Am 2 4m 2 



2 I ' 



where p is the momentum operator and 

1 



(A 2 - m 2 ) 2 



, B = -A, C 



m 2 - 2A 2 
(A 2 -m 2 ) 2: 



(13) 



(ii) 1 in t-channel, 



e -m q r _ g-Ar 



+ 



(m 2 - A 2 )e" Ar 
2A 
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1 



m 2( e ~ mq r _ e -Ar 



+ 



4mi7«2 
1 



U 1 



A 2 )Ae 



-Ar 



4m i7?T,2 



771i771 2 



e -m,r _ e -Ar X 
P 2 + 



777,2- A 2 )e- A ^ 2 

P 



7^(A 2 -m 2 ) 2 

47T777 X 7772 

+ (C-BA 2 )e 



777 i7772 2A 

A(m 9 r + l) e - m9r + 5(Ar + l)e 
pr 



-Ar 



(m) + in s-channel, 



^ 2 (A 2 



\/(r) 



777 



2\2 



4777 17772 



47rr 3 ' 

1 



(14) 



6'p 2 



(a' 2 - 77i 2 )(a 2 - A 2 ) 2 (a' 2 - 777 2 ) 2 (a 2 - A 2 ) s 



26p 2 



{b 2 - 26)p 2 



(a' 2 - m 2 )(a 2 - A 2 ) 3 2a 2 (a' 2 - m 2 )(a 2 - A 2 ) 2 



(15) 



(if) + in u-channel, 
V(r) = 



167rr777 i 777 2 
dr 



re 



-dr 



- dr 



+ 



+ 



cos cr 



2d(c 2 + d 2 ) (c 2 + d 2 ) 2 (c 2 + d 2 ) 



re 



+ 



r 2 e -dr 



+ 



re 



-dr 



-dr 



cos cr 



8d 3 (c 2 + d 2 ) 8d 2 (c 2 + d 2 ) 2d(c 2 + d 2 ) 2 (c 2 + d 2 f (c 2 + d?f 
r sin(cr) 



26 M p 2 



2 cos(cr) 



re 



-dr 



-dr 



+ 



2c{c 2 + d 2 ) 2 {c 2 + d 2 f 4d{c 2 + d 2 ) 2 {c 2 + d 2 ) 



re 



-dr 



- dr 



+ 



+ 



cos(cr) 



2a 2 P 



(16) 



_2d(c 2 + d 2 ) (c 2 + d 2 ) 2 (c 2 + d 2 ) 2 

where the last potential is only displayed in position space for brevity and some small 
relativistic contributions to it are neglected. Here 



(777 2 - mi) 2 , b u 



-(?77 S 



777i 



777 i 777 2 



(17) 



and 



2 2 2 

c =a u -m q 



(m 2 - mi) 2 - m 2 , d = h. - a u = h. - (m 2 - mi) . (18) 

Since the mesons are colour singlets, we do not need to introduce a "confinement" po- 
tential as for the quark-antiquark system. Substituting the potentials into the Schrodinger 
equation, we can solve it and obtain the eigenenergy and wavefunction. As noted previ- 
ously, our interest is in molecular states with small binding energy. In such a molecule 
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the kinetic energy of the constituent mesons in the bound state is also small, so that 
relativistic effects can be neglected, and the non-relativistic Schrodinger equation can be 
used. 



2.3 Decay width 




3 



M 



q 



4 



Figure 1: The decay of meson M. Mesons f and 2 are the constituent mesons in the 
molecular state. Mesons 3 and 4 are the decay products. 

The decay occurs via the triangle diagram shown in Fig. 1. The case of t-channel 
meson exchange is displayed. The s-channel is similar. The diagram represents a Bethe- 
Salpeter (B-S) equation where the internal momentum in the loop is integrated out. The 
two intermediate lines labelled "1" and "2" represent the constituent mesons. In a loosely 
bound molecule they are very close to their mass shell. Thus, instead of the complicated 
B-S integration, we treat the two constituent mesons as free ones, and they transit into 
the final states via an inelastic scattering process. 

When decay via meson exchange takes place, the amplitude is approximately propor- 
tional to the wavefunction (or its derivative) at the origin. This is due to an integration 
effect. When we take the B-S approach to calculate the amplitude via a loop containing 
the bound state with an appropriate kernel (the potential we derived above), we imme- 
diately find that after integrating out the internal momentum, only a wavefunction at 
origin ip(0) remains when the constituent mesons are in relative S-wave. From the for- 
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inula |^(0) | 2 = 2fj, < VItfIV' > @. with the reduced mass of the system, it follows 
that ip(0) is an average over the entire r— region, in other words, effects of all ranges in 
r are involved in ^(O). So even though only the wavefunction at the origin matters, this 
does not mean that only small distance effects are important. For constituent mesons in 
P-wave we know that the coupling is proportional to the momentum. After integration 
in the B-S equation, it turns into the derivative of the wavefunction at origin t/>'(0). 

Let us briefly discuss the legitimacy of the treatment adopted in the decay calculation. 
We find that the validity of the approximation used in the B-S approach depends on the 
size of — being small. From the Heisenberg uncertainty principle, Ipl ~ — . For a 
molecular state with small binding energy E^, the r.m.s. radius r rms is relatively large 
(Eb ~ l/(2/ir 2 m s )). Thus j^- is small as expected. 

The decay can be realized via t- and s-channels and their contributions must be 
summed up. For 0~0~ — > CTCT scattering, exchanging 

(i) + in t-channel, 



gf 



'A 2 — m 



2\ 2 



q 2 — m 2 \ A 2 — q 2 



(ii) 1 in t-channel, 



(19) 



'A 2 -m 2 \ 2 



A 2 - q 2 



{Pi+PzY{P2+P^YD llv {q) 



99 



'A 2 — m 



2\ 2 



q 2 — m 2 \ A 2 — q 2 J 
2M 2 — mf — m\ — m\ — m 2 



(m 2 — m 2 )(m 2 — m|) 
m 2 



+ mi J mj + I 2 + m 2 ym^ + I 



[in 



+ in s-channel, 



(20) 



gf(q — m — iY q m q ) ( A — m 



2\ 2 



(q 2 — m 2 ) 2 + T 2 m 2 \ A 2 — q 2 



(21) 
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where 7713 and are the masses of the outgoing mesons in Fig. 1. Besides these three 
amplitudes, there are those corresponding to CTCT — > 0~1~. We omit their explicit 
expressions. 

The decay width is 

]T Mi\ 2 for 0+ decay, (22) 

i=t,s,u 

Mi\ 2 for 1_ dec ay, (23) 

i=t,s,u 

where % is summed over the s-, t- and u-channel contributions. / = (m 4 + m\ + m\ — 
2m 2 m\ — 2m 2 m\ — 2m 2 m 2 ) 1 ^ 2 / (2m) is the 3-momentum of the outgoing mesons in the 
CM frame. 

3 Numerical results 

Numerical results are obtained by including the a as an exchanged scalar meson, and the 
ou, p and as exchanged vector mesons. There is no low-lying isovector scalar meson, 
and the exchange of the low-lying pseudoscalar mesons tt, i] and rj' is forbidden by parity. 
Vector meson exchange in s-channel is not allowed by conservation of charge conjugation. 
For the ^(1400) and 7T!(1600) systems we also incorporated / (980) and a (980) exchange. 

3.1 / (980) and a (980) 

First we consider / (980) and a (980) which are taken as bound states of KK. Since all 
other parameters are fixed by experimental data and SU(3) flavour symmetry, only the 
off-shellness parameter A is a free parameter. We do not intend to determine its value 
by fitting the mass spectra of / (980) and a (980). In a reasonable range of 1 to 10 GeV, 
we choose a value of A which gives reasonable predictions. We have 

A = 3.4 GeV. 



167rm 3 

wm 2 i 



(ii) iv 



487T7B 5 
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The resultant masses are 



M /o(980) = 10 1 3 MeV M ao(980) = 996 MeV without zero - energy, 
M /o(980) = 999 MeV M ao(980) = 982 MeV with zero - energy, 



and should be compared to the experimental data [16 



M /o(980) = 980 ± 10 MeV M ao(980) = 984.8 ± 1.4 MeV. 

The calculated values are fairly near to the experimental values, given that A is the only 
free parameter we have for the KK system. There is always the possibility to add an 
additional constant term to the potentials that we derived, called the "zero-energy" . The 
term ultimately arises from renormalization when non-perturbative QCD is matched to 
a meson exchange model. In our final mass estimate we included a zero-energy so that 
the calculated values are in closer agreement with experiment. 

The central experimental value of M ao ( 980 ) is a bit larger than that of Mj ( 980 ), but 
our result is inverted. Our results may be changed as follows. The a (980) receives 
no contribution from s-channel a exchange due to isospin conservation. However, the 
/o(980) does. In the case of the / (980), the a intermediate meson constructively interferes 
with t-channel exchange. However, obtaining coupling constants from widths, as we do, 
means that we cannot fix the sign of the amplitude. So there is the possibility that the 
amplitude has an additional —1 between the two channels. This may well change the 
results. Moreover, if we choose different values of A for / (980) and eto(980) and fine-tune 
other parameters, we can substantially improve our fit of the spectra. Since our aim is 
to obtain a general picture with the least parameters, we shall not pursue these further 
refinements. Because parameter changes allow considerable change in binding energy, we 
do not view our calculated binding energy as accurate. Our concern is with an overall 
understanding of the spectra and lifetimes, in order to draw qualitative conclusions. The 
details and subtle deviations are not important. 
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For both / (980) and a (980) 

|^(0)| 2 = 1.2 x 10~ 3 GeV 3 , 

showing that the states are loosely bound. Substituting this value into Eq. [22|, we obtain 

r(/o(980)) — > 7T7r) = 108 MeV s- and t-channels interfere constructively, (24) 
r(/o(980)) — > 7T7r) = 86.5 MeV s- and t-channels interfere destructively, (25) 
r(a (980) -> nri) = 36.3 MeV. (26) 

We display values corresponding to the possibilities of the s- and t-channels interfering 
without or with an additional sign between the two channels. 

The data indicate that the modes /o(980) — > txtx and a (980) — > nr) are dominant, i.e., 
these partial widths can approximately be taken as the total widths. The experimental 
data are 

r(/ (980)) = 40 ~ 100 MeV and r(a (980)) = 50 ~ 100 MeV. 

The values we calculated were not fitted to this data, so that the agreement with experi- 
ment is noticible. 

3.2 7n(1400) and 7^(1600) 

We now calculate the spectra and the decay widths of the 1 h exotic states, assuming them 
to be molecular states. As discussed above, we postulate that 7Ti(1400) and 7Ti(1600) (1 _+ ) 
are molecular states of 7rr](1295) and 7r?7(1440). With the potential derived in Section 2, 
we obtain the following results. 
For tti(1400) 

A = 1.17 GeV,E b = -78 MeV, M^ (1295) = 1354 MeV, |V'(0)| 2 = 3.95 x 10~ 5 GeV 5 . 
For 7Ti(1600) 

A = 1.35 GeV,E b = -14 MeV, M^ {im) = 1521 ~ 1591 MeV, |^'(0)| 2 = 4.18xl0~ 5 GeV 5 
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By contrast to the case for / (980) and a (980), there are several decay modes for 
the exotic states. We can only evaluate a few exclusive modes and later by a general 
argument, we can obtain the order of magnitude of their lifetimes. Below we list several 
partial widths: 



r(7ri(i4oo) 


-> TIT]) 


= 5.28 x 10" 


" 2 MeV, 


r(7ri(i4oo) - 


-> 717]') 


= 6.80 x 10" 


" 2 MeV, 


r(7ri(i4oo) -> 


n~p°) 


= 2.86 x 10" 


" 2 MeV, 


r(7ri(i600) 


-> 7T?7) 


= 4.14 x 10" 


" 2 MeV, 


r(7n(i60o) - 


-> To/) 


= 6.56 x 10" 


" 2 MeV, 


r(7ri(i600) -> 


7T-p°) 


= 3.70 x 10" 


" 2 MeV, 



(27) 

where the calculation was performed for ^1(^-^(1295)) (1400) and 7i"i (^-,7(1440)) (1600). 

In our model, we can estimate the relative decay widths of 7i"i(1600) to 7r6i,7r/i and 
up. The difference is due to the effective couplings and the form factor parameter A. 
With isospin symmetry, we have 

7r&i : 7r/i : i\p ~ 2.6 : 2 : 1. 

One hence obtains various branching ratios for different decay modes, but one can conclude 
that with the hypothesis that the 1 h exotic states are molecular states, tit], 7tt/, Tip, ubi 
and 7r/i are of the same order of magnitude. 

4 Discussion and conclusion 

In this work we consider molecules of two constituent mesons with meson exchange. We 



argued in section |2J] that the picture becomes more accurate if the binding energy of 
the constituent mesons is small compared to their masses. We also argued that low mass 
exchanges dominate in a meson exchange model if the binding energy is small compared 



15 



to the exchange meson mass. In section [2.2| we argued that small binding energy enables 



the use of a non-relativistic Schrodinger equation. In section [2.3| we argued that the 
decay formalism is valid when the binding energy is small with respect to the mass of the 
exchange meson. To summarize, our formalism should be valid when the binding energy 
is small relative to the masses of the constituent and exchange meson masses. The lowest 
mass constituent meson considered in this work is the tt, and the lowest mass meson 
exchanged is found to be the a. 

We investigate the possibility that the observed 1~ + exotic states 7Ti(1400) and 7Ti(1600) 
are molecular states. 

Our strategy is to use the potential derived from field theory to calculate the spectra 
and partial decay widths and see if they can fit the measured data. We derive the 
potential by calculating the elastic scattering amplitude between the constituent mesons 
and then carrying out a Fourier transform to convert the amplitude into the potential in 
configuration space. In our calculations, we obtain the effective couplings from data by 
assuming SU(3) flavour symmetry. Only the free parameter A, which models the off-shell 
effects at the effective vertices, needs to be determined. The inelastic scattering amplitude 
in momentum space is used to calculate the lifetimes. 

To check the plausibility and validity of the approach, we study /b(980) and a (980) 
which are commonly considered as molecular states of KK . The mass spectra and life- 
times of /o(980) and ao(980) are in qualitative agreement with data. Thus we are con- 
vinced that the potential approach is applicable for evaluating the approximate mass 
spectrum and lifetimes of molecular states. 

With the same scheme, we calculate the decay widths of the observed 1 h states 
7r 1 (1400) and 71"! (1600). We assume them to be molecular states of 7r?7(1295) and 7r?7(1440) 
respectively. By fitting the spectra of 7r!(1400) and 7r 1 (1600), we fix the A parameter and 
obtain |?//(0)| 2 which plays a crucial role in the decay width calculation. 

The value of |^'(0)| 2 , and hence the calculated decay widths, are 2 ~ 3 orders of 
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magnitude smaller than needed for consistency with the data. This is the main new 
observation of this work. To understand this we estimate the ratio which governs the 
width of a P-wave molecule relative to that of an S-wave molecule (Eqs. 



R - \M0W ' (28) 

where the subscripts P and S stand for the P and S-waves respectively and m is the mass 
of the molecular state (whose existence can be understood from dimensional analysis). 
As a guide, we note that for a Coulomb potential 

^(0) 

where /i is the reduced mass of the system and a is the Coulomb coupling. For the 
"heavy-light" 1 _+ molecular states consisting of one heavy meson and one light meson 
assumed in this work, the reduced mass is close to the mass of the light meson n, which 
is much smaller than the mass of 7Ti(1400) or 7Ti(1600) (which is close to the mass of the 
heavy meson). Thus there is a suppression factor (/i/m) 2 < 0.01 in the ratio R. Since the 
decay widths of / (980) and a (980) are proportional to ^(O)! 2 , but the decay widths 
of 1~ + molecular states (if they exist) are proportional to \ip' P (0)\ 2 /m 2 , one can easily see 
that a factor of 10 -3 — 10~ 2 would suppress the 1~ + decays. 

There is a close correspondence between the exchange mesons we used, and those 
used in nucleon-nucleon interaction models, e.g. the phenomenologically successful Bonn 
potential | ]32[| . The exchange mesons are the same except for two cases. Firstly, the r] 
and t]' are not employed in the Bonn potential because their couplings to the nucleon 
are small. (For our molecules, these exchanges are not allowed by parity.) Secondly, the 
/o(980) and <p are not considered in the Bonn potential because their large ss content 
implies Okubo-Zweig-Iizuka (OZI) rule violating couplings with the nucleon. For our 
calculations, these mesons have to be considered because the KK molecule has OZI rule 
allowed couplings to ss exchange mesons. The 7r?7(1295) and 7r?](1440) molecules can also 
experience OZI rule allowed / (980) exchange, because the / (980) is not thought to be an 
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almost pure ss state. Unlike the Bonn potential, we do not consider higher order meson 
exchange in this first orientation. 

In our derivations and numerical computations, there are some uncertainties from both 
experimental and theoretical sides. For example, the mass of a— meson ranges from 400 
to 1200 MeV and its decay width is also uncertain. The model has obvious flaws due to 
introduction of phenomenological parameters and ignoring the processes at quark-level. 
Therefore we can never expect very accurate predictions within this framework. If an 
additional factor of 2 ~ 3 exists in the results, it is not surprising. However, since all 
derivations are based on the quantum field theory in a well motivated approximation, 
and parameters are fixed by fitting data, order of magnitude estimates should be valid. 
The only conclusion we can expect is whether the present data can tolerate the molecular 
state scenario for 7i"i(1400) and 7i"i(1600) within experimental errors. 

Since the dissociation rates of the 1~ + molecular states calculated in this work are 
two to three orders of magnitude smaller than the data, the lifetimes of the molecular 
states are fully determined by the lifetimes of the constituent mesons (it, r/(1295) and 
77(1440) resonances). They are not comparable to the measured values of the 7r!(1400) 
and 7Ti(1600) lifetimes. So we conclude that the possibility of 7Ti(1400) and 7Ti(1600) 
respectively being loosely bound 7T77(1295) or 7rr/(1440) molecular states is ruled out by 
the data. 

The decay formalism shows that the dissociation rate of a heavy-light molecule with a 
binding energy much smaller than the masses of the constituent mesons becomes smaller 
the higher the relative angular momentum between the constituent mesons is. The explicit 
calculation suggests that the dissociation rate is already negligible for a P-wave molecule. 
We hence expect the lifetime of a molecule with one or more unit of relative angular 
momentum to be fully determined by the lifetimes of its constituent mesons. Contrary to 
the usual expectation for molecules , such a molecule will have negligible "fall-apart" 
decay by dissociation to the two mesons that it is composed of. Instead, it will decay via 
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the decay of its constituent mesons. A constituent meson will emit at least two particles 
in its decay implying that the decay of the molecule is dominated by its decay to three or 
more particles. This is an unusual decay pattern. 

If an experimental state is found with a lifetime equal to the lifetimes of its con- 
stituents, a probable explanation is that it is a molecule which has a sufficiently low 
binding energy so that decay due to dissociation is small. If the constituent mesons need 
to be in P-wave to obtain the desired quantum numbers of the state, one would usu- 
ally expect a lower-lying S-wave state with different quantum numbers and an increased 
dissociation rate. This could make the S-wave state so broad that it would be difficult 
to observe. We hence predict the possibility of longer living orbitally excited molecular 
states. 
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A Appendix: General formalism 

For interactions among mesons, the early theory is the linear a model and the SU(2) 
Lagrangian is 



L = \[{d,Sf - 2m 2 S 2 ] + (^) 2 i(djf + - XvS* - ±S* + (29) 



where 
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and and a are respectively the pion and a-meson fields. 

Taking into account the interaction with vector mesons, one replaces the differenti- 
ation symbols with covariant differentiation forms. The effective chiral Lagrangian for 



pseudoscalar and vector mesons is j35 

. F 2 



^TriV^UWtf) + -mjjTr(V M V), 



(30) 



where the vector meson fields are 



V„ 



AT 
1 r fi 



V2 



V 



V2 ^ V2 



A* 



P 



Er. _|_ !£lL 1^*0 



(31) 



with A a (a = 1, ...,8) the SU(3) generators. The SU(3) differentiation forms are 



(32) 



where 1^ = v fl + a M and r M = v M — a M are linear combinations of the external vector and 



axial vector a M fields. Also 



U = exp(iA a $ a ) 



(33) 



with $ a the fields of the pseudoscalar meson octet. 

From recent studies, the a— meson may be an independent field, so that we retain the 
interaction in Eq. In our final expression, we keep the leading a— related terms and 
the leading terms in the pseudoscalar-vector interactions, while we ignore all axial vectors 
because they are much heavier. We also extend the Lagrangian to include the iso-scalars 
7/(1295) and 7/(1440). The much simplified Lagrangian reads 

g*PP + g'v i P^P k + ..., 



J eff 



(34) 



where k guarantees a proper SU(2) combination, and the ellipsis denotes the smaller 
contributions which are neglected in our calculations. We note that the symmetry reduces 
from SU(3) to SU(2). 
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